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Abstract:Here, the model of non-equilibrium model with two states (— 1,+1) 
and a noise q on simple square lattices proposed for M.J. Oliveira (1992) following 
the conjecture of up-down symmetry of Grinstein and colleagues (1985) is stud- 
ied and generalized. This model is well-known, today, as Majority- Vote Model. 
They showed, through Monte Carlo simulations, that their obtained results fall 
into the universality class of the equilibrium Ising model on a square lattice. In 
this work, we generalize the Majority- Vote Model for a version with three states, 
now including the zero state, (—1,0,-1-1) in two dimensions. Using Monte Carlo 
simulations, we showed that our model falls into the universality class of the spin- 
1 (—1,0,-1-1) and spin-1/2 Ising model and also agree with Majority- Vote Model 
proposed for M.J. Oliveira (1992) . The exponents ratio obtained for our model 
was 7/z/ = 1.77(3), ji/v = 0.121(5), and 1/z/ = 1.03(5). The critical noise obtained 
and the fourth-order cumulant were qc = 0.106(5) and U* = 0.62(3). 
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1 Introduction 

The Ising model PQ [2] has been used during long time as a "toy model" for 
diverses objectives, as to test and to improve new algorithms and methods of 
high precision for calculation of critical exponents in Equilibrium Statistical 
Mechanics using the Monte Carlo method as Metropolis [3J, Swendsen-Wang 
[1], Wang-Landau [S] algorithms. Single histogram [B] and Broad histogram 
[7| methods. The Ising model was already applied decades ago to explain how 
city populations segregate [13j, how a school of fish aligns into one direction 
for swimming ^ or how workers decide whether or not to go on strike |9l 
[To] . In the Latane model of Social Impact [llj the Ising model has been 
used to give a consensus, a fragmentation into many different opinions, or a 
leadership effect when a few people change the opinion of lots of others. To 
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some extent the voter model of Liggett [12] is an Ising-type model: opinions 
follow the majority of the neighbourhood. All these cited model and others 
can be found out in 

Realistic economics models of tax evasion appear to be necessary, be- 
cause tax evasion remains to be a major predicament facing governments 
[151 [m [m [IS]- Experimental evidence provided by Gachter [T2j indeed 
suggests that tax payers tend to condition their decision regarding whether 
to pay taxes or not on the tax evasion decision of the members of their 
group. Frey and Torgler [20] also provide empirical evidence on the rele- 
vance of conditional cooperation for tax morale. Following the same context, 
recently, Zaklan et al. [211 E2] developed an economics model to study the 
problem of tax evasion dynamics using the Ising model through Monte-Carlo 
simulations with the Glauber and heatbath algorithms (that obey detailed- 
balance equihbrium) to study the proposed model. G. Grinstein et al. [23] 
have argued that nonequilibrium stochastic spin systems on regular square 
lattices (SL) with up-down symmetry fall into the universality class of the 
equilibrium Ising model [22] ■ This conjecture was confirmed for various 
Archimedean lattices and in several models that do not obey detailed bal- 
ance [211 1251 1261 123 1211 129] . The majority- vote model with two states (MV2) 
is a nonequilibrium model proposed by M.J. Oliveira in 1992 and defined 
by stochastic dynamics with local rules and with up-down symmetry on a 
regular lattice shows a second-order phase transition with critical exponents 
/3, 7, V which characterize the system in the vicinity of the phase transition 
identical [23 [30] with those of the equilibrium Ising model [1] for regular 
lattices. Lima et al. [31] studied MV2 on Voronoi-Delaunay random lattices 
with periodic boundary conditions. These lattices posses natural quenched 
disorder in their connections. They showed that presence of quenched con- 
nectivity disorder is enough to alter the exponents f5/v and 7/z/ from the 
pure model and therefore that quenched disorder is a relevant term to such 
non-equilibrium phase-transition which disagree with the arguments of G. 
Grinstein et al. [23] . 

Recently, simulations on both undirected and directed scale-free networks 
[321 [221 [2S1 [Ml [33 [Ml [Sg, random graphs [IHl [H] and social networks 
[121 [131 [H], have attracted interest of researchers from various areas. These 
complex networks have been studied extensively by Lima et al. in the context 
of magnetism (MV2, Ising, and Potts model) [13 [10 |13 [IS [13 [50] , econo- 
physics models [221 [M] and sociophysics model [52] ■ Lima [53] make an anal- 
ysis of tax evasion dynamics with the Zaklan model on two-dimensional SL 
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using MV2 for their temporal evolution under different enforcement regimes. 
They showed that the MV2 model also is capable to control the different levels 
of the tax evasion as it was made by Zaklan et al. [22] using Ising model on 
various structures: SL, Voronoi-Delaunay random lattice, Barabasi-Albert 
(AB) network and Erdos-Renyi (ER) graphs; we discuss the resulting tax 
evasion dynamics. Perhaps, this is the first application of MV2 model to 
a real system, in this case applied to the economy. The generalization to 
a three-state majority-vote model(MV3) on a regular SL was considered by 
[5H [55] . where the authors found Qc = 0.117(1). The resulting critical ex- 
ponents for this non-equilibrium MVS model are in agreement with the ones 
for the equilibrium three-states Potts model [SB], supporting the conjecture 
G. Grinstein et al. [23]. Recently, Melo et al. [SZ] also studied the MVS 
model[5ll |55], now, on ER graphs [IQ]. In MVS model on ER graphs they 
found that the critical noise Qc is a function of the mean connectivity z of the 
graph. The critical exponents ratio 7/z/, /3/z/ and l/u was calculated for var- 
ious values of connectivity. In this, Melo et al.[57] suggest a future work on 
the two- and three-state Potts model on random graphs would be of interest 
in order to provide a direct comparison with our results in light of the con- 
jecture by Grinstein et al. [23], which states that reversible and irreversible 
models with same symmetry belong to the same universality class. 

In the present work, we propose a non-equilibrium model with three states 
(—1, 0, +1), called the three-state majority-vote model (MVS) in two dimen- 
sions as a generalization of MV2. Our main goal is to check the hypothesis of 
Grinstein et al. [23] described above. The remainder of our paper is organised 
as follows. In section 2, we present our MVS model and detail of the Monte 
Carlo simulation and calculations used in the evolution of physical quantities 
of MVS. In section S we make an analysis of simulations performed in the 
previous section and discuss the results obtained. And finally in section 4, we 
present our conclusions from the results obtained of MVS model and present 
several perspectives work to be done at present. 
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2 Model and simulation 



On a SL where each site of the lattice is inhabited, at a time step, we consider 
the MV3 defined by a set of "voters" or spin variables a taking the values ±1 
and 0, situated on every node of the SL with N — sites. The evolution is 
governed by single spin-fiip like dynamics with a probability Wi of i-th spin 
flip given by 



Wi{a) 



h-{l-2q)aiS(f2'^i+s) 



(1) 



where Six) is the sign ±1 of x if a; 7^ 0, S[x) = if a; = and sum runs over 
the number z = 4 oi nearest neighbors of i-th spin. The control parameter 
< q < 1 plays the role of the temperature in equilibrium systems and 
measures the probability of parallel aligning to the majority of neighbors. 
It means, that given spin i adopts the majority sign of its neighbors with 
probability q and the minority sign with probability {1 — q) . 

To study the critical behavior of the model we define the variable m = 
mi — 777,3 if 777 1, 7772, and 7778 are the numbers of sites with ctj = +1, and 
— 1, respectively, all normalized by A^. In particular, we are interested in the 
magnetization M, susceptibility x the reduced fourth-order cumulant U 

M{q) ^ {\m\), (2) 

X(g) = N ((777^) - (771)^) , (3) 

U{q) ^ 1 - (m^)/ {?>{m^?) , (4) 

where (■ ■ ■) stands for a thermodynamics average. The results are averaged 
over the Nmn independent simulations. These quantities are functions of the 
noise parameter q and obey the finite-size scaling relations 

M = L-^'^fm{x), (5) 

X = L^/Vx(^), (6) 

f = LVV.(.), (7) 

where j3, and 7 are the usual critical exponents, fi{x) are the finite size 
scaling functions with 

x^{q- q,)L'/^ (8) 
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being the scaling variable. Therefore, from the size dependence of M and x 
we obtained the exponents (i/u and 7/z/, respectively. The maximum value 
of susceptibility also scales as U^^'^ . Moreover, the value of q* for which x 
has a maximum is expected to scale with the system size as 

q* = q, + bL-'/" (9) 

where b ^ 1. Therefore, the relations may be used to get the exponent 
l/u. We performed Monte Carlo simulation on the SL with various systems 
of size L = 8, 16, 32, 64, 128 and 256. It takes 1 x 10^ Monte Carlo steps 
(MCS) to make the system reach the steady state, and then the time averages 
are estimated over the next 2 x 10^ MCS. One MCS is accomplished after 
all the N spins are investigated whether they flip or not. We carried out 
Nrun = 1000 to 10000 independent simulation runs for each lattice and for a 
given set of parameters {q,N). 

3 Results and discussion 

In Fig.l, we show the dependence of the magnetization M on the noise param- 
eter q, obtained from simulations on SL with lattice size L = 8, 16, 32, 64, 128 
and 256 with [L x L = N) sites. The shape of M{q) curve, for a given value 
of N, suggests the presence of a second-order phase transition in the system. 
The phase transition occurs at the value of the critical noise parameter q^. 

In Fig. 2 the corresponding behavior of the susceptibility x is presented. In 
Fig. 3 we plot the Binder's fourth-order cumulant U for different values of the 
system size N. The critical noise parameter qc is estimated as the point where 
the curves for different system sizes intercept each other [59]. The critical 
noise obtained was qc = 0.106(5) and "universal charge" U* = 0.62(3). In 
Fig. 4 we plot the dependence of the magnetization M* = M{qc) versus the 
linear system size L. The slopes of curves correspond to the exponent ratio 
/3/z/ according to relation (5). The obtained exponent is /3/z/ = 0.121(5). The 
exponents ratio 7/1/ are obtained from the slopes of the straight lines with 
7/1/ = 1.77(3) for SL, as presented in Fig. 5 and obtained from the relation 
(6). To obtain the critical exponent l/i/, we used the scaling relation (9). 
The calculated values of the exponents 1/z/ = 1.03(5) . See Fig. 6. 

The results of simulations together with data for ER graphs and AB 
networks with mean connectivity z = 4 are collected in Tab. 1. 
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Figure 1: The magnetization M versus the noise parameter g, for L = 8, 16, 
32, 64, 128 and 256 size lattice for SL. 

4 Conclusion 

We presented a very simple non-equilibrium MVS on SL. On this lattice, the 
MVS shows a second-order phase transition. Our Monte Carlo simulations 
demonstrate that our MVS model agree with the universality hypothesis of 
Grinstein et al. |2Sj . 

Finally, we remark that the critical exponents 7/1^, P/v and 1/v for MVS 
on SL presented in the Tab. 1 are similar the exponents from the = 1/2 
and S = 1 Ising model 2D [H ED], obtained of analytical results, and also 
with the results obtained, through Monte Carlo simulation, by M. J. Oliveira 
[27] and Kwak et al. [61]. Therefore, our results agree with the conjecture by 
Grinstein et al.[23], which states that reversible and irreversible models with 
same symmetry belong to the same universality class. Future works using 
MVS model are already in progress on undirected and directed Small-world 
networks and on ER random graphs. 
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figure 1. 



Susceptibility versus q for SL. Tlie same simulation design as in 
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